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Abstract. 

We discuss the physics of backreaction-driven accelerated expansion. Using the exact 
equations for the behaviour of averages in dust universes, we explain how large- 
scale smoothness does not imply that the effect of inhomogeneity and anisotropy 
on the expansion rate is small. We demonstrate with an analytical toy model how 
gravitational collapse can lead to acceleration. We find that the conjecture of the 
accelerated expansion being due to structure formation is in agreement with the general 
observational picture of structures in the universe, and more quantitative work is 
needed to make a detailed comparison. 
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1. Introduction 



Evidence for acceleration. There is a large body of observational evidence supporting 
the claim that the expansion of the universe has accelerated in the recent past, and 
may be accelerating today. This conclusion has been bolstered by the verification of the 
prediction of the location of the baryon acoustic peak in the matter power spectrum [1], 
in a convincing demonstration of concordance. In addition, the worrisome feature that 
nearby and distant populations of type la supernovae used to have different absolute 
magnitudes and were both individually consistent with deceleration has disappeared 
with new and better data [2] (though see [3,4]). 

The ACDM model, where the acceleration is driven by vacuum energy (or the 
cosmological constant, which is the equivalent modification of gravity) agrees well with 
most observations, with the notable exception of the low CMB multipoles [5-9] (it has 
also been argued that cluster observations support a non-accelerating universe [10]). 
However, given the lack of theoretical understanding about the parameters of the ACDM 
model (notably the vacuum energy density), it is a phenomenological fit rather than a 
well-founded theory, and its success does not rule out the possibility that quite a different 
model can also be a good fit to the data. (The values obtained for the parameters of a 
cosmological model by fitting to observations should not be mistaken for measurements, 
as model selection studies show; see e.g. [3,4,11].) In particular, while the observation 
that there is accelerating expansion seems robust, the nature of the acceleration is not 
well constrained. In the ACDM model, the transition to acceleration is gradual, but a 
rapid transition is not ruled out [12,13]. In fact, from the SNIa data it is difficult to say 
anything beyond that the universe has accelerated in the recent past, even whether the 
expansion is still accelerating [3,4,14]. 

Keeping to the assumption that the universe is completely homogeneous and 
isotropic, any explanation of the acceleration has to involve either a medium with 
negative pressure or modified gravity. Such models in general, and the ACDM model 
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in particular, suffer from the coincidence problem: why does the acceleration happen 
around a redshift of unity, at around 10 billion years? In other words, why are we seeing 
a very particular phase in the evolution of the universe, when the inferred energy density 
of the source driving the acceleration has recently become equal to the energy density of 
matter? The clearest qualitative change in the late-time universe is the formation of non- 
linear structures. It therefore seems a natural possibility that the observed deviation 
from the prediction of homogeneous and isotropic cosmological models with normal 
matter and gravity could be related to the known breakdown of the assumption that 
the universe is homogeneous and isotropic (rather than to a speculated failure in the 
description of the matter content or the theory of gravity). 

The inhomogeneous universe. One possible avenue is trying to explain the observations 
without having any accelerated expansion. Cosmological information is borne to us by 
light along null geodesies (apart from information carried by neutrinos and cosmic rays). 
The standard analysis of light propagation assumes that the universe is perturbatively 
near a homogeneous and isotropic Friedmann-Robertson- Walker (FRW) model, which 
is manifestly not true on scales smaller than 70-100 h~^Mpc [15-17]. 

It is therefore possible that the propagation of light would be affected by 
the inhomogeneities and/or anisotropics in a way that looks like acceleration when 
interpreted in the context of an FRW model. Studies of the Lemaitre-Tolman-Bondi 
(LTB) model [18] (see [19] for a review), the spherically symmetric dust solution of 
the Einstein equation, have demonstrated that the effect of inhomogeneity on the 
luminosity distance can mimic acceleration [20,21] (see [21] for more references). Even 
though spherical symmetry is a questionable assumption for the entire universe, it could 
be a good first approximation for the local region. In any case, one would expect 
a qualitatively similar effect to be present also in more realistic and less symmetric 
spacetimes [22] - arguably, the effect of dumpiness could even be stronger when there is 
less symmetry. The effect of inhomogeneity and anisotropy on the luminosity distance 
has also been studied in perturbed FRW models [23,24]. 

An explanation of the apparent acceleration in terms of inhomogeneity and/or 
anisotropy could solve the coincidence problem, since inhomogeneity and anisotropy 
become important only in the late-time universe. However, inhomogeneity and 
anisotropy affect different observations in different ways, and it would require an odd 
coincidence for all the various indicators of expansion rate (SNIa luminosity distances, 
the cosmological microwave background (CMB) anisotropics, large scale structure (LSS), 
and so on) to be affected in a way that would be consistently interpreted as acceleration 
when fitting to a FRW model. 

One proposed possibility is that we live in an underdense region, a 'Hubble bubble' 
[25] (for more references, see [21]). In this proposal, the local matter density today 
is Qmo ~ 0.15-0.35, as indicated by local observations [26], while the global value is 
Qm = 1- The SNIa luminosity distances as well as the difference between the local and 
global values of the expansion rate could be explained in terms of inhomogeneity, while 
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a global model with no acceleration can fit most other observations, including the CMB 
and LSS [27] (though it is not supported by studies of the local and global expansion 
rate [28]). However, it would be difficult to explain the baryon acoustic peak [29]: one 
would have to appeal to inhomogeneities (or features in the primordial power spectrum) 
to supply a pattern that by coincidence happens to fit the expectations of an accelerating 
FRW model. 

One way to phrase the issue is that cosmological observations involve a larger 
number of a priori independent parameters than the ACDM model. Therefore the 
ACDM model implies relations between observationally independent parameters. (For 
an early discussion of cosmological observations in an inhomogeneous and anisotropic 
spacetime which makes the issue transparent, see [30].) It is not surprising that models 
with more degrees of freedom, such as the LTB model which involves two arbitrary 
functions, could fit the data as well as FRW models. However, it would be an unlikely 
coincidence for them to also produce the same relations between observables as FRW 
models. This is generally true for any models where the explanation of the luminosity 
distances is decoupled from the explanation of the low matter density, baryon acoustic 
peak and so on (such as mixing of photons with axions [31] or with the gauge bosons of 
a new f/(l) gauge group [32]). 

The fitting problem. While the FRW scale factor has been very successful in fitting 
observations, it is difficult to understand the matter content implied by the FRW 
equations which relate the scale factor to the energy-momentum tensor. Note that 
the fact that the mean properties of the universe are well described by an overall scale 
factor does not imply the stronger statement that the scale factor evolves according to 
the FRW equations, since the universe is not completely homogeneous and isotropic. 

The idea that the average behaviour of inhomogeneous and/or anisotropic 
spacetimes is in general different from the behaviour of homogeneous and isotropic 
spacetimes goes back to at least 1963 [33]. The first comprehensive discussion was given 
in 1983 by George Ellis [34], who called the task of finding the smooth metric which 
best fits the real clumpy universe the fitting problem. The influence of inhomogeneity 
and/or anisotropy on the average behaviour is also known as backreaction [23,35-82]; 
see [19,54] for further references, in particular early ones, and [83] for an overview. 

The idea that perturbations with wavelengths smaller than the Hubble radius 
could lead to acceleration for the scale factor (as opposed to merely mimicking the 
appearance of acceleration via changing null geodesies) was studied in the context 
of linear perturbation theory in [54] (the possibility had been earlier touched upon 
in [44,48]; see also [50]). The calculation was then extended to second order [23,58], 
and it was suggested that linear perturbations with wavelengths much larger than the 
Hubble radius could lead to acceleration [23,58-60]. It is now agreed that this is not 
possible [61-63, 65, 69j3. As for perturbations smaller than the Hubble radius, there 

I Even though super-Hubble perturbations do not contribute to acceleration in a dust universe, 
they could lead to deceleration during inflation driven by a scalar field or a cosmological constant, 
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is no acceleration to at least second order in perturbation theory [54,58,69,77,78]. 
If inhomogeneity and anisotropy are to explain the observed acceleration, the only 
possibility is via non-linear sub-Hubble perturbations, that is, the process of structure 
formation, as proposed in [50,54]. 

It has been analytically shown in the LTB toy model how backreaction of non-linear 
perturbations can modify the Hubble law [57], and acceleration has also been numerically 
demonstrated in the LTB model [76,79,80], but the physical meaning of inhomogeneity- 
and anisotropy-driven acceleration and the connection to structure formation has been 
unclear. We will discuss the relation between homogeneity and isotropy, the overall scale 
factor, the FRW metric and the equations which describe the average expansion of the 
universe. We will then look at an exact toy model of structure formation, explain the 
physics of acceleration driven by inhomogeneity and anisotropy, and note that structure 
formation involves a preferred time near the observed acceleration era. In particular, we 
will clarify two apparent paradoxes of backreaction-driven acceleration: how the average 
expansion of a manifold can accelerate even though the local expansion rate decelerates 
everywhere, and how collapse implies acceleration. These issues were earlier discussed 
in the brief essay [81]. 

In section [2] we discuss the assumptions underlying FRW models, and go through 
the derivation of the Buchert equations which describe the average behaviour of an 
inhomogeneous and/or anisotropic dust spacetime. In section [3] we consider an exact toy 
model where gravitational collapse produces acceleration, discuss how this mechanism 
may operate in the real universe, and compare this picture with some observations and 
simulations of structures in the universe. In section H] we summarise the situation with 
regard to the conjecture that the observed acceleration is due to backreaction. 

2. Smoothness and variance 

2.1. The Friedmann-Robertson-Walker assumptions 

The three assumptions. The assumption of homogeneity and isotropy that underlies 
the Friedmann-Robertson- Walker models of cosmology can be broken down into three 
distinct parts. 

Assumption 1: FRW scale factor. The observables characterising the mean properties 
of the universe can be computed from an overall scale factor. 

Assumption 2: FRW dynamics. The overall scale factor evolves according to the FRW 
equations. 

Assumption 3: FRW + perturbations. Deviations from homogeneity and isotropy evolve 
according to linear perturbation theory around the average. 

and supcr-Hubblc scalar field perturbations left over from inflation could still be important today 
[37, 39, 49, 52, 53, 67, 68, 72, 73, 75]. 
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Since the universe is not exactly homogeneous and isotropic, it is clear that the above 
assumptions do not hold exactly. However, there is usually no attempt to quantify 
the deviation, and the first two assumptions are often conflated. A notable exception 
is the program of observational cosmology (and related work) by George Ellis and 
collaborators, which aimed at formulating cosmological theory in a manner that does 
not involve a priori assumptions about homogeneity and isotropy and that is as close to 
the observations as possible [84-87]; see [88,89] for an overview. 

We will look at the above assumptions from the slightly different point of view 
of backreaction studies, where the emphasis is not so much on establishing the level 
of inhomogeneity and isotropy of the universe as on quantifying their effect on the 
expansion rate and other observables. Let us discuss the three assumptions in turn. 

Assumption 1: FRW scale factor. Observations indicate that the universe is 
statistically homogeneous and isotropic on large scales. It has recently become possible 
to directly establish from the Sloan Digital Sky Survey, for the first time, the average 
homogeneity of the universe by looking at the fractal dimension of the point set of 
luminous red galaxies [15,16]. The related homogeneity scale has been quantified as 
70-100 /;,~^Mpc, though analysis based on the morphology of structures indicates a 
homogeneity scale that is larger than 100-200 /i~^Mpc [17]. Dividing the observational 
volume into 10 regions with individual volume 2x10^ (/i~^Mpc)'^ (corresponding to 
a ball with radius ~ 170 /i~^Mpc), the density variance is 7% in the redshift range 
0.2 < z < 0.35, quantifying the degree of statistical homogeneity in present observations. 
The largest structure known is the Sloan Great Wall at 420 /z~^Mpc, which has 
superseded the old 240 /i~^Mpc Great Wall. These sizes are 14% and 8%, respectively, of 
the Hubble radius; in the Einstein-de Sitter universe (the spatially flat matter-dominated 
FRW model) this would be 7% and 4%, respectively, of the visual horizon. Structures 
this large are rare, and the typical size of observed collapsing structures and voids is 
^ 20-40 /i~^Mpc, of the order 10~^ to 10~^ of the visual horizon [90-92] (though some 
simulations suggest that a significant fraction would be larger [93] or smaller [94]). 

Statistical isotropy is in turn supported by the high degree of isotropy in the CMB, 
along with the 'almost Ehlers-Geren-Sachs theorem', which states that a universe where 
the CMB looks almost isotropic everywhere is almost FRW on large scales [95, 96] 
(see [97,98] for discussion and caveats). The applicability of the 'almost EGS theorem' 
to the real universe is, in fact, somewhat unclear. The proof of the theorem requires the 
assumption that the expansion rate is positive everywhere, which is not valid in the real 
universe unless scales where structure formation by gravitational collapse is relevant are 
smoothed over. Even then, it is unclear how the strict limits (essentially given by the 
CMB anisotropy of 10^^) on spatial variation of the exansion rate and other observables 
can be reconciled with the observed (and theoretically expected) differences of order 
one in the expansion rate and energy density in non-linear regions (we will discuss the 
observations in section [331) • Nevertheless, it seems intuitively clear that the isotropy of 
the CMB (coupled with the assumption that we do not occupy a preferred position in 
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the universe) indicates a high degree of average isotropy in the geometry on large scales. 

Statistical homogeneity and isotropy show that a description in terms of an overall 
scale factor could be consistent. In physical terms, because the size of typical non-linear 
regions is small, light rays coming to us cosmological distances pass through several such 
regions on the way to us, and the differences could be expected to average out. However, 
the small size of non-linear regions does not prove that a description in terms of an overall 
scale factor is necessarily correct. We receive almost all cosmological information along 
null geodesies, and from the fact that inhomogeneities and anisotropies are small when 
averaged over large scales it does not follow that a description in terms of an overall scale 
factor will correctly capture the physics of light propagation, as discussed in section [H 

The FRW metric is conformally fiat, so the Weyl tensor which embodies the non- 
local effects of gravity vanishes, while the Ricci tensor, determined by the local matter 
distribution, is non-zero. However, photons in the real universe only occasionally 
encounter matter and mostly travel in vacuum, where the Weyl tensor is non-zero but the 
Ricci tensor vanishes (neglecting the small contribution of the microwave and neutrino 
backgrounds). The FRW description is exactly the opposite of the real situation, and 
it is not obvious why it would correctly capture the physics of the passage of light in 
the real inhomogeneous and anisotropic universe. From a theoretical point of view, this 
issue has not been satisfactorily settled [99-102]; for an overview, see [83] and for further 
references, see [21,83]. 

A related issue is that cosmological models are usually discussed in terms of 
hypersurfaces of proper time, whereas observations are mostly made along the past 
lightcone. This issue was considered in detail in the program of observational cosmology 
[85]. In FRW models this distinction is not important, and they enjoy (for monotonous 
expansion) a simple one-to-one correspondence between time and redshift. In an 
inhomogeneous and/or anisotropic universe, the issue is more complicated, and is related 
to the problem of choosing the hypersurface of proper time [49,53,54,57]. Looking only 
at an average scale factor has been criticised in [71] (see also [74]) on the grounds that 
one obtains different behaviour for different choices of time slicing. 

In FRW models with only a single fluid, there is a preferred time coordinate given 
by the proper time measured by observers comoving with the fluid. This notion can 
be straightforwardly extended to inhomogeneous and/or anisotropic dust spacetimes 
which do not have any symmetries. If the matter consists of irrotational dust, then such 
hypersurfaces of constant proper time are everywhere orthogonal to the fluid flow lines, 
and the situation is analogous to the FRW case. If vorticity is present, the hypersurfaces 
of proper time will not mesh together to fill the spacetime [103,104], and the situation 
is more complicated [63,65,96]. (During inflation driven by a cosmological constant or 
a scalar field, the issue of time slicing and observables is more involved [39,49,53,68].) 
In any case, averaging does involve a loss of covariance; see [105] for discussion. 

The approximation that the matter is an irrotational, pressureless ideal fluid will 
necessarily break down at small scales once structure formation has started [106,107], as 
otherwise collapsing structures could not stabilise, given that the local expansion rate 
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of irrotational dust can never increase. In other words, structure formation and the 
associated vorticity need to be considered and may be important for backreaction. As 
we will see, zero vorticity is an important mathematical assumption for the definition 
of averages and an overall scale factor. However, from a physical point view one would 
not expect the vorticity associated with structure formation to make a quantitative 
difference for the expansion of the universe. The issue should be carefully considered, 
but we shall simply assume that the small-scale breakdown of the picture of matter as 
an irrotational ideal fluid with zero pressure is not important. 

Note that the problem of time and the averaging hypersurface is not an artifact 
of inhomogeneous and/or anisotropic models, but a feature of the general relatistic 
description of the real universe. It is a virtue of the observational cosmology and 
backreaction approaches that they make these issues explicit and provide tools for 
studying and quantifying them, unlike FRW models. 

A seemingly worrisome aspect of the hypersurfaces of constant proper time is that 
observables averaged over them depend on regions outside the visual horizon, since 
the hypersurfaces extend beyond the past lightcone. It may seem unphysical that an 
observational quantity would change if we changed the definition of the hypersurface in 
regions beyond our past lightcone. However, if the universe is statistically homogeneous 
and isotropic on large scales, we do not have the freedom to independently adjust the 
hypersurface of proper time in the regions inside and outside our past lightcone, since 
the state of matter and geometry inside and outside the lightcone is required to be 
statistically identical for the same proper time. (See [84,86] for discussion of homogeneity 
in cosmology.) This puzzle is therefore just the homogeneity and isotropy problem, 
which is solved (or at least alleviated [108]) by inflation. The statistical equality of 
widely separated regions was set up in the early universe when they were in causal 
contact, so the particle horizon is much larger than the visual horizon. (For clarification 
of the different horizons, see [109, 110].) In spacetimes which are not statistically 
homogeneous and isotropic, the averaging procedure would not be expected to be useful, 
and a description in terms of an overall scale factor would probably not make sense. 
For example, this would be the case if the size of typical non-linear structures were a 
sizeable fraction of the visual horizon. 

We will simply extend the FRW notions of proper time and scale factor in the most 
straightforward manner (as we will discuss in section 12. 2p and assume that redshift is 
related to the scale factor in the same way as in FRW models. Dependence on proper 
time could then be expressed in terms of redshift, as in the FRW case, so that the 
problem discussed above is not apparent. Of course, this is simply a matter of rewriting 
the equations in a manner that makes the assumption of statistical homogeneity and 
isotropy less transparent. 

Properly, one should derive the quantitive conditions under which the scale factor 
approximation is valid and see to which extent they are realised in the universe. 
This would involve defining redshifts, luminosity distances and other observables using 
null geodesies in the inhomogeneous and anisotropic universe, and determining under 
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which conditions these reduce to the quantities defined with the overall scale factor. 
As discussed above, the issue of average light propagation in inhomogeneous and/or 
anisotropic spacetimes has not been satisfactorily settled from the theoretical point of 
view. 

Nevertheless, from an observational point of view, the scale factor ansatz has 
been very successful in fitting a range of observations. This success is all the more 
remarkable as different observations depend on null geodesies in a different manner. For 
example, there are several different definitions of the expansion rate (such as the volume 
expansion rate, the rate of deviation between neighbouring geodesies, the expansion rate 
inferred from the luminosity distance and so on), which agree for a homogeneous and 
isotropic space but differ when inhomogeneities and/or anisotropies are present. Yet, 
expansion rates inferred from different observations agree to within ^ 20% [28]. The 
approximation of considering only the scale factor has worked very well in practice, 
whatever its theoretical status, and the observed statistical homogeneity and isotropy 
on large scales at least guarantees that this approximation is consistent. 

Note that even small corrections to the approximation of looking only at the scale 
factor can be interesting for the light they shed on the average behaviour. For example, 
the dipole of the angular power spectrum of the inhomogeneous luminosity distance 
in a linearly perturbed FRW model yields a direct measure of the Hubble parameter 
as a function of redshift [24]. Also, even if the passage of light through cosmological 
distances, and thus many non-linear regions, is well described in terms of a scale factor, 
this does not rule out the possibility that local structures could affect null geodesies in 
a way that is not captured by that approximation. In particular, the anomalies of the 
low CMB multipoles [5-9] could be related to the local breakdown of the description in 
terms of a linearly perturbed FRW metric [111,112]. 

(Another distinct issue, which we will not discuss, is the "dressing" of cosmological 
parameters, i.e. the feature that the usual interpretation of observations does not 
account for the geometrical inhomogeneities and/or anisotropies in the averaging domain 
when considering observables [46,51,55].) 

Even if the effects of inhomogeneity and anisotropy cancel on large scales so that 
the FRW scale factor is a good average description, this does not mean that the FRW 
metric would be a good approximation of the average geometry. The reason is that the 
FRW metric also contains the spatial curvature, which is assumed to be homogeneous 
and isotropic. In a general spacetime, the spatial curvature is inhomogeneous and 
anisotropic, and does not evolve like the FRW spatial curvature, even on average. The 
reason is that the evolution of non-linear regions with positive spatial curvature and 
those with negative spatial curvature is different, and they are not correlated so as to 
produce the FRW behaviour. In particular, even if a universe is perturbatively close 
to spatial flatness early on, this condition is not necessarily preserved once density 
perturbations become non-linear. 

In mathematical terms, this is related to the fact that the FRW evolution oc 
of the spatial curvature in terms of the scale factor a arises from the integrability 
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condition between the Raychaudhuri equation and the Hamiltonian constraint. The 
two receive different corrections from the inhomogeneities and anisotropics, so the 
general integrabihty condition does not agree with the FRW case and the average spatial 
curvature does not in general evolve like a~^. From a physical point of view, regions with 
negative spatial curvature expand faster than regions with positive spatial curvature, so 
one would expect that they will come to dominate the volume and the average spatial 
curvature will become negative. The non-FRW evolution of the spatial curvature and 
the competition between overdense and underdense regions is central to the proposed 
backreaction explanation for accelerated expansion, and we will discuss these issues in 
detail in sections 12.21 and 13. 1[ The behaviour of spatial curvature is at the heart of 
the distinction between the scale factor being a good description and the scale factor 
following the FRW equations, an issue to which we now turn. 

Assumption 2: FRW dynamics. As discussed above, we will simply assume that the 
average properties of the universe can be described in terms of an overall scale factor, 
without deriving the conditions under which this assumption is valid. In contrast, when 
considering the dynamics, we will write down the exact equations governing the evolution 
of the scale factor, quantify the domain of validity of the FRW equations and discuss 
the impact of the corrections. Indeed, while there is strong observational support for 
the FRW scale factor, the time evolution given by the FRW equations is quite different 
from what is observed, unless the equations are amended by adding a source term with 
negative pressure or by modifying gravity. 

It has been suggested that in order to emphasise the difference between the 
geometry and the dynamics, the names Robertson-Walker would be associated with 
the assumption that the geometry is approximately homogeneous and isotropic, while 
the stronger assumption that the dynamics of the scale factor is given by the Einstein 
equation applied to a completely homogeneous and isotropic metric would bear the 
names Friedmann-Lemaitre [88, 89]. (As discussed above, the present situation is 
slightly different in that we are not even assuming that the geometry is described by a 
homogeneous and isotropic metric, simply that we can use an overall scale factor.) 

Simply inserting the overall scale factor into the Einstein equation is not the 
correct way to find the dynamical equations which it satisfies. Instead, one should 
insert the full inhomogeneous and/or anisotropic metric into the Einstein equation and 
only then take an average, since the average behaviour of an inhomogeneous and/or 
anisotropic spacetime is not the same as the behaviour of the corresponding smooth 
spacetime (where "corresponding" means that the smooth and average quantities have 
the same initial conditions). In other words, the average properties of an inhomogeneous 
and/or anisotropic spacetime do not satisfy the Einstein equation. The fact that taking 
an average metric and plugging it into the Einstein equation and plugging the real 
metric into the equation and then averaging are not equivalent is sometimes expressed 
by saying that because the Einstein tensor G^^, is non-linear in the metric, we have 
{G^uiOa/s)) 7^ G^^{{ga/3)), whcrc () stands for averaging. It would be more accurate 
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to say that the problem arises because time evolution and averaging do not commute. 
(The statement about averaging the Einstein equation is anyway rather heuristic because 
tensors cannot be straightforwardly averaged on curved manifolds; though see [113].) 
This is the origin of the fitting problem discussed in section [1] 

The standard assumption is that even after non-linear structures have started 
forming, the average evolves according to the FRW equations if smoothing on the scale 
of the non-linearity is performed. However, the equations for the mean expansion which 
are actually derived, rather than assumed, do not bear out this expectation [36,41,47]. 
As we will discuss in section \2.2\ inhomogeneities and anisotropics affect the behaviour 
of the scale factor, and neither statistical homogeneity and isotropy nor the small size of 
the inhomogeneous and anisotropic regions (relative to the visual horizon) is a sufficient 
condition for recovering the FRW behaviour. 

A different reason why the FRW equations could be invalid even though the universe 
is very homogeneous and isotropic on large scales has also been advanced. It has been 
suggested that the breakdown of the approximation of treating the matter as an ideal 
fluid involves negative pressure which might explain the acceleration [50]. An argument 
against the treatment of matter as an ideal fluid goes as follows. Small-scale processes 
in the universe are in general not thermodynamically irreversible but instead produce 
entropy. Since there is no production of negative entropy, this entropy generation 
does not vanish upon averaging. In contrast, in the approximation of an adiabatically 
expanding ideal fluid, the entropy of the universe is constant, and such small-scale effects 
are completely absent. In numerical terms, a single 3 x 10^ solar mass black hole, such as 
the one at the center of our galaxy, has an entropy of 10^°, of the same order of magnitude 
or more as the total entropy ascribed to the observable universe in the ideal fluid picture. 
Supermassive black holes are abundant in the universe, so the entropy associated with 
them alone completely overwhelms the entropy associated with the adiabatic fluid. On 
the other hand, in the average description of an ideal fluid spacetime it may look as if 
the entropy was increasing, even though there is no local entropy production [56]. The 
issue of gravitational entropy is tied up with coarse-graining and thus backreaction, and 
is not well understood [83] . Whether these problems of the ideal fluid description are 
important for the expansion rate is not clear. 

A related issue is whether one can neglect the influence of the structure within 
stabilised regions on the overall expansion, i.e. whether one can continuously and 
consistently "renormalise" the scale of the stable regions which are treated as particles 
of the dust fluid [114]. 

We will not study these two issues, and keep to the assumption that deviations 
from the dust behaviour are not important for the overall expansion. 

Assumption 3: FRW + perturbations. The region of validity of the assumption that 
inhomogeneities and anisotropics on a FRW background evolve according to linear 
perturbation theory is well-known, and its breakdown at the end of the linear regime 
is well understood. As long as the density contrast 6 = (p — (p))/ (p) of a perturbation 
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is small, it evolves according to linear perturbation theory. In a spatially flat matter- 
dominated FRW background, the growing mode is proportional to the scale factor, 
5 (X a. As 6 becomes of order ±1, the linear approximation breaks down. An 
overdensity will collapse and the density contrast will grow faster than in the linear 
regime, whereas an underdensity will grow more slowly than in the linear regime, 
asymptotically approaching emptiness. 

For intermediate regime perturbations which have not yet gone non-linear but are 
nevertheless small enough to be located inside a non-linear structure, one would expect 
the evolution to depend on the environment and not just on the overall average expansion 
of the universe. For example, the evolution of a 20 h~^Mpc radius perturbation would 
be expected to be different inside a 50 h~^Mpc void or a 200 h~^Mpc wall. Such 
perturbations, which are larger than the size of typical structures and yet fit inside 
non-linear regions, are by definition untypical, and one can argue that their effect on 
the mean evolution of perturbations is small. 

However, it is not obvious that linear perturbation theory around the average in a 
space which is highly inhomogeneous and anisotropic correctly captures the evolution 
even for perturbations with wavelength longer than the size of the largest structures. 
Regions with large overdensities or underdensities presumably contribute differently to 
the evolution of the long-wavelength modes, and it is not clear that these effects would 
cancel to give the same answer as linear perturbation theory around the average. 

The separation into background and perturbations is more involved in a universe 
with large inhomogeneities and anisotropics than in the FRW case, and in section 
12.21 below we will show that linear perturbations do not in general satisfy the same 
perturbation equation as in FRW models. We will not discuss the issue further, but 
in a realistic model of backreaction, the sensitivity of perturbations to smaller scale 
inhomogeneities and anisotropics should be looked at in more detail, particularly when 
these have a large effect on the background. For work on perturbations in a backreaction 
context, see [38,40,42,43,48]. 

2.2. The Buchert equations 

The metric and the local equations. We assume that the matter content of the universe 
can be described as dust, i.e. a pressureless ideal fluid. We further assume that the 
dust is irrotational (i.e. the vorticity is zero). Then the metric can be written in the 
synchronous gauge [103, 104] 



where t is the proper time measured by observers comoving with the dust and ^^"'g^j^t, x) 
is the metric on the hypersurface of constant t. The Einstein equation reads 



ds^ = -dt^ + '^^^g,Ax'dx^ , 



(1) 



87rGj\fTf^i, 



(2) 
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where G^^(t, x) is the Einstein tensor, Gjy is Newton's constant, T^^{t, x) is the energy- 
momentum tensor, p{t,x) is the dust energy density and = (1,0) is the velocity of 
comoving observers. 

We wish to find the equations for average quantities. Since only scalars can be 
straightforwardly integrated on a curved manifold (though see [113]), we should project 
([2]) to obtain a set of scalar equations. In addition to u'^ and g'^'^, we have the covariant 
derivative available. From these we can build three independent rank two tensors 
to project witli§|, so the Einstein equation ([2]) yields the following three exact, local, 
covariant scalar equations [103,104,115] 

e + \e'^ = - AttGnP - (3) 

V =8vrG^p-l(3)i? + a2 (4) 
p^Qp =0, (5) 

where a dot stands for derivative with respect to t, 9(t,x) = {^/^^^g)^^dt{^/^^^g) is 
the expansion rate of the local volume element, a'^{t,x) = 2a'^^aij > is the scalar 
built from the shear tensor (Tjj, and ^^^R{t,x) is the Ricci scalar of the hypersurface of 
constant t (i.e. the spatial curvature). The acceleration equation ([3]) is known as the 
Raychaudhuri equation, and (jll) is the Hamiltonian constraint. 

The price for reducing the Einstein equation down to a set of scalar equations is 
that the system is not closed: there are three equations for four independent variables. 
Essentially, the propagation of the shear tensor (or equivalent ly, of the Ricci tensor on 
the hypersurface of constant t) does not reduce to a scalar equation. The integrability 
condition between ([3]) and (jlj) reads 

dti'^'^R) + le = 2 dta' + 49a^ , (6) 

so specifying either the shear or the spatial curvature fixes the other. 

Note that no approximations have been made: the equations (I3])-(I5]) are exact for 
irrotational dust, with arbitrarily large density variations. 

Deriving the Buchert equations. We are interested in the evolution of average 
quantities: specifically, we want to know how the average expansion rate behaves. Our 
discussion follows the original derivation by Buchert [41]. When (and only when) the 
vorticity is zero, the rest spaces of constant proper time of comoving observers mesh 
together to form a family of hypersurfaces which fills spacetime. The spatial average of 
a scalar quantity / is then straightforwardly defined on these hypersurfaces as 

jd^x,/i^y^)fit,x) 

where the integral is over the hypersurface of constant t. An important property of the 
averaging ([7]) is that it does not commute with time evolution, 

dt{f) = {dtf) + {fe)-{f){e) . (8) 

§ Since V^Gnv vanishes identically, V^V^ does not give anything new compared to m^V. 
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In order to describe the average evolution and compare it to FRW models, we have 
to define a scale factor. The simplest extension of the notion of an overall scale factor 
to an inhomogeneous and/or anisotropic spacetime is to define it with the volume of the 
hypersurface of constant t: 

where the normalisation has been chosen as a (to) = 1 at some time to (which we shall 
take to be today). In words, a{t)^ is the volume of the hypersurface of constant t (up 
to the usual multiplicative constant). We could equivalently define the scale factor with 
the average of the volume expansion rate, by 3a/a = (6). We will also use the notation 
H = a/a. 

By taking the average ([7]) of the scalar equations ([SJ-Q and commuting the time 
derivatives as shown in we obtain the equations satisfied by the scale factor first 
derived by Thomas Buchert in 1999 [41]: 

3- = -47rG,v(p) + Q (10) 
a 

3^ = 87rG^(p)-lf)i?)-iQ (11) 

d,{p) + 3-{p) = 0, (12) 
a 

and the integrability condition between the average Raychaudhuri equation ( fTOl) and 
the average Hamiltonian constraint ffTTl) . analogous to IQ, reads 

9^((3)^) + 2^((3)/2) = -Q-6^Q, (13) 
a a 

where the backreaction variable Q is a new term compared to the FRW equations, 
containing the effect of inhomogeneity and anisotropy: 

Q^l{{e')-{e)')-2{a') . (14) 

The Buchert equations fll0p - fll3p are exact for the averages when matter consists of 
irrotational dust. (The Newtonian limit was derived by Buchert and Ehlers in 1995 [36] 
and the case with non-zero pressure by Buchert in 2001 [47].) The backreaction variable 
Q consists of two terms: the variance of the expansion rate and the shear. Shear is also 
present in the local equation (j3]) and decelerates expansion. In contrast, the variance 
is only present in the averaged equations, arising from the the non-commutation of 
averaging and taking a time derivative as shown in ([8]), and acts to accelerate the 
expansion rate. The presence of this term makes it possible for the average equations 
to display behaviour which is qualitatively different from the local behaviour. In the 
Newtonian limit, Q can be written in terms of the Minkowski functionals, which are 
a statistical measure of the morphological properties of cosmic structure, relating Q 
directly to structure formation [44]. 
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Recovering the FRW equations. As with the system (E])-®, there are only three 
independent equations in f|T0l) - f|T3l) for the four independent functions a, (p), Q and 
Physically this means that different inhomogeneous and/or anisotropic 
spacetimes can evolve differently even if they have the same average initial conditions. 
While the equations (fT0l) - (fT3l) cannot be solved, they can be used to check whether a 
given scale factor can result from backreaction [65]. 

Specifying one more function (or relation between the four functions) leads to a 
soluble system. In particular, in the limit when the shear and the variance of the 
expansion rate are small compared to the contribution of the energy density, the Buchert 
equations reduce to the FRW equations. The integrabihty condition (fT3ll between (fTOll 
and (fTT]) then leads to the standard behaviour {^^^R) oc a~^. (Since the integrabihty 
condition ties the evolution of spatial curvature to the backreaction variable Q, the 
average spatial curvature evolves like {^^^R) cx only in the FRW limit or in the 
special case when Q oc a~^, in agreement with the discussion of assumption 1 in section 
12.11 ) Conversely, if the shear or the variance of the expansion rate are large compared 
to the contribution of the energy density in a large fraction of space, the FRW equations 
are not a good approximation for the average behaviour. This is the way to derive the 
FRW equations and quantify their domain of validity. 

The average behaviour also reduces to the FRW equations when the shear and the 
variance of the expansion rate cancel, even if they are not small. (For examples of 
such solutions, see [35].) In the Newtonian limit, this happens generically for periodic 
boundary conditions [36] and for spherically symmetric spaces [42]. (With Minkowski 
functionals, the latter result can be understood to follow from the property that the 
backreaction variable Q in the Newtonian limit measures the deviation of morphology 
from that of a ball [44].) This cancellation is also present in relativistic perturbation 
theory when expanding to second order, as shown in [54] and recently rederived in [77]. 
This is not evidence towards a theorem that backreaction would not in general affect the 
acceleration, as claimed in [77] . Instead, it is related to the vanishing of backreaction for 
periodic boundary conditions in the Newtonian limit. It was correctly identified in [64] 
that in the perturbative expansion the Newtonian terms are those with the largest 
number of spatial gradients, as they are accompanied by the largest number of powers 
of the speed of light. Therefore, the terms with the highest number of spatial gradients 
at each order in perturbation theory vanish for periodic boundary conditions (which are 
implicit in the use of Fourier decomposition). Since the highest number of gradients 
at order in perturbation theory is 2A^, the second order term with four gradients 
vanishes upon averaging. However, at fourth order there is no reason for the term with 
six derivatives to vanish, and this term becomes large after non-linear structures start 
forming, so the perturbation expansion is expected to break down. 

An important point is that the expansion described by the Buchert equations 
f|T0|l - f|T3l) does not reduce to the FRW behaviour simply when the spatial size of 
the inhomogeneities and anisotropies is small. If the shear or the variance of the 
expansion rate is comparable to the contribution of the energy density in a sizeable 
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fraction of space, the behaviour will deviate from the FRW case (barring the sort of 
cancellation discussed above), regardless of the size of the individual inhomogeneous 
and/or anisotropic regions. The variance of the expansion rate is non-negative, and 
only vanishes if the expansion rate is completely homogeneous: the contributions from 
different inhomogeneous regions cannot cancel to zero. The same is true for the shear. 
(Of course, the shear and variance from different regions could cancel each other.) The 
Buchert equations f|T0l) - f|T3l) make explicit and quantify the statement made in section 
12.11 that large-scale statistical homogeneity and isotropy does not guarantee that the 
dynamics of the scale factor follows the FRW equations. 

Acceleration without acceleration. The average Raychaudhuri equation (fTO!) shows that 
if the variance of the expansion rate is large enough compared to the shear and the 
energy density, the average expansion accelerates, even though the local expansion 
rate decelerates at every point according to the local Raychaudhuri equation ([3]). The 
possibility of acceleration is a property of the averaged system which is not present in 
the local behaviour, and is due to the non-commutation of averaging and time evolution. 

Since Q contributes positively to the acceleration (ITOl) . but negatively to the Hubble 
rate (fTTl) . it might seem that negative curvature is needed to balance the negative 
contribution of Q to the Hubble rate, as claimed in [65]. In fact, it is possible to have 
acceleration even when the spatial curvature is positive, as all that is needed is that the 
sum —{^^^R) — Q increases (i.e. becomes less negative). This just means that {^^^R) 
has to decrease (i.e. —{^^^R) has to increase) faster than Q is growing. Acceleration 
necessarily involves (non-FRW) spatial curvature, as f|T3|) shows: if {^^^R) (x a~^, we 
have Q oc a~^, and there is no acceleration. This feature demonstrates the difference 
between having a description in terms of a scale factor and the FRW metric being valid, 
discussed in the context of assumption 1 in section 12.11 

One consequence of the fact that local expansion can only decelerate is that once 
a shell of matter has started collapsing, it cannot turn around and stabilise. It follows 
that the formation of stable structures involves vorticity and/or breakdown of the dust 
approximation [106, 107]. The contribution of vorticity to the local acceleration ([3]) 
(which we did not include) is always positive, so it cannot disappear upon averaging, 
as the contribution from each region has the same sign. This is also the case for 
shear, and the two can cancel each other. Indeed, in stabilised regions, the positive 
contribution of vorticity has to exactly balance the negative contributions of the shear 
and the energy density to produce net zero acceleration. One would thus naively expect 
our approximation of neglecting vorticity to lead to a lower bound on the acceleration. 
However, vorticity would also complicate the averaging proceduce, as discussed in section 
[H so the issue should be carefully studied. 

We will discuss the physical meaning of the average acceleration and the relation 
to structure formation in section [3l but let us first complete the overview of the FRW 
assumptions by looking at assumption 3 concerning the behaviour of small perturbations 
in an inhomogeneous and/or anisotropic universe. 
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The evolution of perturbations. We will briefly consider perturbation theory to see what 
are the corrections to the FRW picture. For linearly perturbed FRW models, one writes 
the equations of motion as the sum of the background part plus a small perturbation. 
Discarding terms beyond first order in the perturbation and taking the spatial average 
leads to the equations for the FRW background (since the average of the perturbation is 
taken to vanish). Deducting these from the perturbed equations then gives the evolution 
equation for the perturbations. 

When backreaction is important, the same procedure does not work, as the 
difference between the local equations ©-([S]) and the averages (fT0|) - (fT2l) is large by 
definition. We will do the closest thing, which is to separate the local terms into a 
'large' part which determines the average behaviour and a 'small' part which does not 
contribute to the averages. We write 9{t, x) = 9o(t, x)+A9(t, x), and assume that |^o| ^ 
|A6'|. We similarly split p{t,x) = po{t,x) + Ap{t,x), a'^{t,x) = a^{t,x) + A{a'^){t,x) 
and ^^''R(t,x) = ^^^Ro(t,x) + A{^^^R)(t,x). The assumption that the 'large' part is 
responsible for the average behaviour, (6) = (6*0), implies that {A6) = 0, and similarly 
for the other quantities. We then insert this ansatz into the local equations ([3])-([5]) and 
keep only terms up to linear order in the 'small' quantities. Unlike in the FRW case, 
there is no rigorous way to separate the perturbation from the background, as both 
have spatial dependence. A similar issue arises in deriving the perturbation equations 
for 'tilted' cosmological models, i.e. models where the fluid velocity is not everywhere 
normal to the hypersurface of constant proper time [116]. (If we had non-zero vorticity, 
the model would necessarily be tilted.) If we simply assume that the equations for the 
'large' and 'small' parts decouple, then the 'large' parts satisfy the local equations ([3])- 
([5]) and their average gives the Buchert equations (fT0l) - (fT2l) . For the 'small' terms, we 
then obtain the following evolution equation 



where 6 = Ap/po, along with the consistency condition {9qA9) = 0. In addition to the 
shear source term on the right-hand side, f|T5|) differs from its FRW counterpart in that 
9o,po are position-dependent. 

It is possible to obtain evolution equations where the average expansion rate and 
energy density appear instead of the local quantities. Let us assume that S is separable, 
6{t,x) = D[t)h{x)] in the linearly perturbed FRW case, this would correspond to 
considering pure growing or decaying modes. Dividing (fTS!) by h{x) and averaging, we 
have 



The evolution equation ( fT6ll agrees with the linearly perturbed FRW casq||l, apart from 
the shear source term (note that (p) oc by f|T2l) ). In the spatially fiat FRW 
limit, we recover the standard behaviour D oc a for the growing mode. Let us write 

II For a given a{t); the evolution of a (and thus the evolution of D) as a function of time will in general 
be different, since it is governed by the Buchert equations and not by the FRW equations. 




(15) 




(16) 
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D(t) = a(t)d{t) in order to analyse the deviation from this hmiting case. Using equations 
f fTOj) and f|TT]) we can write f|T6|) as 



If the expansion accelerates, the spatial curvature will become more negative, 
tending to decrease d, just like acceleration and negative spatial curvature act against 
the growth of perturbations in FRW models. The shear term can have either sign: 
writing a^^ = + Act*-' for the shear tensor, we have A((j^) = Aa^Aaij. 

Just as the equations for the averages are not closed, neither are the 
perturbation equations. The perturbation equations require more information about the 
inhomogeneities and anisotropies (specifically about the shear perturbation) in addition 
to that needed for closure of the average equations. Once this information is supplied, 
the backreaction equations for the averages and the perturbations are analogous to 
modified gravity: there is no new energy component, but the relation between the 
average energy density and the average expansion rate, as well as between the averages 
and the perturbations, is different from the FRW case. For the averages, any model of 
modified gravity can of course be formally written as a general relativistic model with 
extra sources, and vice versa. Likewise, in the case of backreaction one can write the 
average behaviour in terms of a scalar field, the "morphon" [82] . It would be interesting 
to see how far the analogy can be extended to the perturbations. 

We will not discuss perturbation theory further. In a realistic model of backreaction 
the issue should be studied in detail, in order to be able to compare with observations of 
CMB and LSS. For work on perturbations in a backreaction context, see [38,40,42,43,48]. 

3. Acceleration from collapse 

3.1. A toy model for backreaction 

Breakdown of the FRW approximation. The proposal [54] that structure formation 
leads to accelerating expansion implies in the context of the Buchert equations ( JTOl) . 
( IT4|) that the relative variance of the expansion rate is of order one. This behaviour 
has been numerically demonstrated in specific examples of LTB models [76,79,80]. 
However, the physical meaning of having average acceleration even though the local 
expansion decelerates everywhere has not been clear [71] (see also [74]). 

We will clarify the physical meaning of the acceleration with a simple toy model, 
and show that it is not unreasonable for structure formation to involve a variance in the 
expansion rate that is large enough to produce acceleration. A shorter treatment was 
presented in the essay [81], and the role of collapsing regions has also been discussed 



As mentioned in section 12.11 when discussing assumption 3, the breakdown of 
perturbation theory at 5 ~ ±1 is well understood. The simplest model used to 
describe the non-linear evolution is the spherical collapse model for overdense regions 





in [79]. 
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(see e.g. [117,118]), and the equivalent for underdense regions (the appendix of [119] 
has a useful summary of both cases). The model consists of a spherically symmetric 
density perturbation embedded in a FRW universe (surrounded by an empty region 
to make the mean density agree with the background), studied in the Newtonian limit. 
The perturbation behaves on average like an independent FRW universe with positive or 
negative spatial curvature in the case of an overdense or underdense region, respectively. 
In terms of the Buchert equations f|T0|) - f|T3|l . this comes about because the backreaction 
variable Q vanishes for spherical symmetry in the Newtonian limit [42] . (For some work 
on the general relativistic case, see [57,80].) 

According to the spherical collapse model, the expansion of an overdense domain 
slows down with respect to the background until the structure stops expanding, turns 
around and starts collapsing. The spherical collapse model indicates that matter 
collapses to a singularity in a finite time, whereas in practice structures stabilise at 
a finite radius, usually taken to be half the radius at turnaround. The spherical collapse 
model is not a good description of the final stages of collapse, as departures from 
spherical symmetry are important for real structure formation, and collapse amplifies 
asymmetries. There are more accurate and realistic treatments of collapse, see e.g. 
[43,45,106,120], but as we are mostly interested in qualitative features, the spherical 
collapse model will be an adequate description of non-linear density perturbations. 

What is not generally appreciated is that just as linear perturbation theory around 
the FRW universe breaks down as perturbations become non-linear, the FRW universe 
itself breaks down as a description of the average behaviour when the non-linear 
perturbations occupy a sizeable fraction of space, as discussed in section 12.21 We want 
to look at this breakdown of FRW equations with a simple model of structure formation, 
in analogy with the description of the breakdown of perturbation theory in the spherical 
collapse model. 

Two-region toy model of structure formation. In the real universe, structure formation 
consists of small overdensities and underdensities developing into stable structures with 
fixed density and voids which are constantly becoming emptier, respectively. We will 
consider the simplest possible toy model of structure formation, with two separate 
spherically symmetric dust regions, one overdense and one underdense. We will consider 
the Newtonian limit, so that the regions evolve according to the spherical collapse model. 
The regions are taken to be disjoint, and we ignore their embedding into the whole 
space. We are then essentially just comparing two FRW universes. The same kind of 
toy model was qualitatively discussed in [71]; we will take a more detailed look and add 
an understanding of the physics behind the equations. 

We denote the scale factors of the regions by oi, 02 and the corresponding Hubble 
parameters by Hi, H2, where region 1 is underdense and region 2 is overdense, so 
Hi > H2. Since the regions are disjoint, the total volume is simply the sum of the 
volumes af and a^, and the overall scale factor is a = (af -|- a^)^^^. The overall Hubble 
and deceleration parameters are (they can be computed from (fTOl) . ffTTj) . f|T^ or directly 
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from the definition of a) 

3 3 

H = -^^H^ + --^^H2 = H^{l-v + vh) (18) 

^ ^-v vh^ v{l-v){l-hf 

^~ H^a'^^^l-v + hvY^'^^l-v + hvY {l-v + hvf ' ^ ^ 
where gi, q2 are the deceleration parameters of regions 1 and 2, and we have denoted the 
fraction of space in the overdense region by v = a\/ {a\ + a\) and the relative expansion 
rate of the two regions hy h = H2/H1. 

The Hubble rate f|T8l) is simply the volume-weighted average of the Hubble rates in 
regions 1 and 2. Not so for the deceleration parameter: in addition to the first two terms 
in (fT9|) . there is a third term related to the variance of the expansion rate, corresponding 
to the backreaction variable Q in the average Raychaudhuri equation (fTU]) . While qi, q2 
are positive or zero, the last term is always negative, corresponding to the fact that Q 
is positive. 

For simplicity, we take the underdense region to be completely empty, so ai oc t. 
Region 2 behaves like a closed FRW universe, with 02 oc 1 — cos0, t (x (p — sin 0, where 
the parameter (j) is called the development angle. The overdense region starts expanding 
from the big bang singularity at = and slows down until it turns around at = vr and 
starts collapsing, finally shrinking to zero size and infinite density at = 27r. In practice, 
overdense regions stabilise at fixed size and density. In the spherical collapse model, this 
is often implemented by hand at = 37r/2, when the radius of the structure is half the 
radius at turnaround. We will therefore follow the evolution only until = 37r/2. 

There is one free parameter in this toy model, the relative size of the two regions 
at some given time. Taking this time to be at turnaround at = vr and denoting the 
fractions of space in regions 1 and 2 at that moment by /i = 1 — /2, /2, we have 

/2vr^(l — cos0)^ 



8(1 - /2)(0 - sin 0)3 + /27r3(l - cos , 
sin 0(0 — sin 0) 



(20) 



(1 — COS0)^ 

Inserting fl20|) into f|T9|) . it is easy to establish that the acceleration can be positive. 

In figured] (a) and (b), we have plotted the deceleration parameter q and the Hubble 
rate multiplied by t, as dimensionless measures of the acceleration and the expansion 
rate, respectively. In addition to the toy model, we show the behaviour in the ACDM 
model, just to make the qualitative features of backreaction-driven acceleration easier to 
grasp by comparison; the toy model is not meant to be taken seriously in a quantitative 
sense. We have chosen the value of the free parameter to be /2 = 0.3, so that the value 
of g at = 37r/2 in the toy model approximately equals the value in the ACDM model 
when Qm = 0.3, Q\ = 0.7. In figure [H (c) we show the density parameters for matter, 
spatial curvature and the backreaction variable Q (defined in section [331 after fl2T|) . fl22l) ) 
in the toy model. Note that negative corresponds to positive spatial curvature, and 
vice versa. 



Accelerated expansion from structure formation 



20 




12 3 



(a) (b) (c) 

Figure 1. The evolution of the toy model as a function of the development angle </). 
(a): The deceleration parameter q in the toy model (blue, solid) and in the ACDM 
model (red, dash-dot), (b): The Hubble parameter multiplied by time, Ht, in the 
toy model (blue, solid) and in the ACDM model (red, dash-dot), (c): The density 
parameters of matter, spatial curvature and the backreaction variable Q in the toy 
model. Red (solid) is n„i, blue (dash-dot) is fi^ and green (dash) is fig. 

Figure [1] (a) shows that the expansion accelerates, particularly after region 2 starts 
collapsing at = tt. It may seem paradoxical that gravitational collapse induces 
acceleration. However, the explanation is simple. Initially, the expansion is similar 
to the Einstein-de Sitter case, with q = 1/2, Ht = 2/3. The contribution of the 
overdense region VH2 gradually slows down the average expansion rate H. The relative 
volume V occupied by the overdense region decreases monotonously, and eventually it 
becomes so small that the contribution of the underdense region begins to dominate 
and the expansion rate, as measured by Ht, grows. The related change from positive 
to negative spatial curvature is clearly seen in figure [1] (c). This effect is particularly 
pronounced and easy to understand once the overdense region has started collapsing: 
then the Hubble rate H2 is actually negative, and its contribution shrinks rapidly as 02 
contracts, so the average expansion rate rises. (The fraction of volume in the overdense 
region at = 37r/2 is only v ~ 0.01.) The increase in the absolute value of the collapse 
rate cannot compensate for the decrease in volume, as the FRW Hubble law shows: 
afHi = al^/SnGNpiQai^ /3 — Kia^'^ . This also means that the results are not sensitive 
to the diverging behaviour in the final stages. The contribution of the collapsing region 
decreases as it approaches the singularity, and the mean quantities would remain finite 
even if we followed the system to the end of the collapse. 

Note that we have first averaged over the internal behaviour of each region, and 
only then taken the average over the two regions. This neglects the variance within 
each region, thus underestimating the total variance by the volume-weighted sum of the 
individual variances. Since the regions are spherically symmetric, this contribution 
to the backreaction variable Q (defined in ([HI)) is exactly cancelled by the shear 
contribution, so the two-step method does give the right Q. 
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General lessons. The reason that the average expansion can accelerate even though 
the local expansion decelerates everywhere is that the growth of the relative volume 
occupied by the faster expanding regions contributes to the average acceleration. This 
makes the physical content of the Buchert equation (fTOj) clear: the larger is the variance 
of the expansion rate, the wider is the difference between the fastest and the slowest 
expanding regions, so the more rapidly the relative volume of the fastest region can 
grow, and the stronger is the acceleration (assuming that the shear is not too large). 

It is also transparent why acceleration from backreaction involves decreasing spatial 
curvature, as indicated by (fTOj) . (ITT!) . (fT4|) and shown by figure [1] (c): the more 
negatively curved a region is, the faster it expands, and the larger its relative volume 
will become. Therefore the most negatively curved regions will come to dominate the 
average curvature. For underdense FRW regions. Hi is above 2/3 and approaches unity 
from below. So, the decrease of spatial curvature is expected to be accompanied by Hi 
approaching 1, as seen in figured] (b). 

While the rise of Hi beyond the Einstein-de Sitter value of 2/3 is due to the 
underdense region, it is the overdense region which is essential for acceleration, i.e. the 
fall of q. (It is possible, but harder, to have acceleration without a collapsing region.) 
Unlike in the ACDM model, these two effects are distinct. For example, if we replaced 
the empty void obeying ai (x t with the Einstein-de Sitter universe for which ai oc 
there can still be acceleration because of the collapse, even though Ht is at most 2/3. 
The important thing is that Ht first slows down, so that it can later rise rapidly due to 
the collapsing region. 

In fact, a more slowly expanding underdense region can even lead to stronger 
acceleration. For example, if instead of ai oc t we were to take ai oc which IS 
arguably closer to the behaviour of real voids [121], q would be more negative than in 
the case when region 2 is empty or an Einstein-de Sitter universe. The reason is that 
there are two competing effects: if the expansion of the faster expanding region is too 
slow, the contrast with the overdense region won't be strong enough to give a large 
variance, while if it's too rapid, the faster expanding region will dominate the Hubble 
rate before the overdense region can have any impact. 

The fact that gravity is attractive (for matter satisfying the strong energy condition) 
implies that the local expansion rate is bounded from above, as can be shown by 
integrating the Raychaudhuri equation (fTOl) as an inequality. The same holds for 
the mean expansion rate, resulting in the bound Ht < 1 for acceleration driven by 
backreaction [65]. In contrast, the collapse rate is not bounded from below, so collapsing 
regions can lead to arbitrarily rapid change in the expansion rate, and q can become 
arbitrarily negative. In fact, q can diverge to negative infinity in a finite time. In 
the two-region toy model, for sufficiently large /2 the negative contribution VH2 will at 
some point equal the positive contribution [1 — v)Hi, so H passes zero from above and 
q diverges to positive infinity. When H later passes zero from below on its way back to 
positive values, q will diverge to negative infinity. 

Such behaviour is in contrast to FRW models, where Ht can grow without limit, but 
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g < — 1 requires violating the null energy condition (or the modified gravity equivalent). 
Acceleration with q < —1 is consistent with observations, and at present there is no 
statistically significant evidence either for it or against it [3,4,12]. 

3.2. Backreaction in the real universe 

Start of the acceleration. The toy model studied above shows that gravitational collapse 
is intimately associated with accelerating expansion, and makes transparent the physical 
content of the Buchert equations f|T0l) - f|T3l) . Obviously, the real universe does not consist 
of two disjoint spherically symmetric regions. Let us clarify which features of the toy 
model are expected to be relevant for the real universe. The density distribution of the 
universe is characterised by a hierarchy of perturbations nested within perturbations. 
In the past the perturbations were small on all scales, so the universe was well described 
by linear perturbation theory around the FRW universe. (And it has been shown that a 
solution of the linearised equations indeed corresponds to a nearly FRW solution of the 
the full equations [122].) In the simplest models of infiation, the primordial spectrum of 
density perturbations is adiabatic and nearly scale-invariant (at horizon entry), and for 
typical models of supersymmetric weakly interacting dark matter there is a cut-off in 
the power spectrum at small scales due to collisional damping and free-streaming [123]. 

The density perturbations grow logarithmically during the radiation-dominated era 
and linearly during the matter-dominated era. For typical supersymmetric dark matter, 
the first generation of perturbations becomes non-linear and forms bound structures 
around a redshift of 40-80 [123]. For other viable dark matter candidates such as axions 
[124], light dark matter [125] or right-handed neutrinos [126], structure formation may 
begin at a different time due to differences in the details of the perturbation spectrum, 
but the qualitative picture is expected to remain the same. For a comprehensive analysis 
of dark matter candidates, see [127]. 

The structures which collapse and stabilise (or form voids) are part of larger scale 
perturbations which are in turn undergoing the process of slowing down and collapsing 
(or speeding up and becoming emptier), and so on, with ever larger non-linear overdense 
structures and underdense voids forming over time. Since total mass is conserved, the 
formation of high-density regions is accompanied by the formation of larger underdense 
regions. 

One possibility is that the gravitational collapse associated with the formation of the 
first generation of structures would lead to acceleration. Since acceleration smoothens 
inhomogeneities and anisotropies and impedes structure formation, backreaction-driven 
acceleration will eventually end. This is implicit in the limit Ht < 1, which rules 
out eternal acceleration [65]. However, in the case of the real universe the physical 
reason for the end of the acceleration is not as clear as in the two-region toy model, 
though according to the perturbation equation (fT7|) the negative spatial curvature 
will work against structure formation. After acceleration ends, inhomogeneities and 
anisotropies can become important again, as perturbations are nested inside each other. 
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with new modes constantly entering the horizon. This could lead to oscillations between 
deceleration and acceleration. Such oscillations are not ruled out observationally, and 
might be detectable in the next generation of experiments [128]. Oscillation between 
acceleration and deceleration is a possible solution to the coincidence problem, as one 
could observe acceleration in the recent past at all times after the start of structure 
formation. 

The question of whether acceleration does start with the collapse of the first 
generation of structures and then undergo oscillations should be addressed in a realistic 
model. At first sight, such early acceleration seems unlikely because backreaction 
involves dynamical spatial curvature, as we have discussed. It seems likely that 
significant spatial curvature from redshift ~ 40-80 onwards would have been apparent 
in observations, particularly in the CMB [129]. However, precisely because the spatial 
curvature does not behave in the same way as in FRW models, the CMB bounds cannot 
be straighthforwardly applied; see section 13.31 for discussion of spatial curvature and 
observations. 

Oscillations aside, backreaction offers another way of addressing the coincidence 
problem. In addition to the time when the first generation of objects forms, structure 
formation involves a second preferred time, near the era ~ 10 billion years when 
acceleration has been observed. 

The process of hierarchical structure formation is not entirely self-similar even 
when the primordial spectrum of perturbations is scale- invariant. In addition to the 
cut-off scale due to collisional damping and free-streaming, there is at least one other 
scale present, related to the transition from radiation domination to matter domination 
(the physics of dark matter or inflation can of course involve further scales). The 
scale keq corresponding to the wavenumber of the mode which entered the horizon at 
matter-radiation equality is imprinted on the perturbation spectrum because metric 
perturbations which entered the horizon during the radiation-dominated era {k > keq) 
are damped relative to those which enter during matter domination {k < k^q). In 
the real universe, the matter-radiation equality scale is k~g ^ 60-160 h'^Mpc (for 
0.15 < Qrno < 0.35 [26] and the updated value Hq = 62.3 ± 5.2 km/s/Mpc [28]). 
Therefore the amplitude of the metric perturbations grows with increasing wavelength, 
asymptotically approaching the value A a; 10~^ set by inflation (or some other process 
in the primordial universe). 

There is a factor k'^/{aH)'^ in going from the metric perturbations to density 
perturbations, and the suppression for small wavelengths is approximately l/k"^ up to a 
logarithm, so there is only a logarithmic difference in the amplitude of the density 
perturbations at small wavelengths, as mentioned above. As a result, the size of 
structures which are just starting to collapse relative to the Hubble radius, Rni/{aH)~^, 
grows rapidly at first and rises monotonically, slowing down as structure formation 
proceeds and saturating at the value ~ \/A once all perturbations which entered the 
horizon during the radiation-dominated era have collapsed. (Here Rni{t) is the scale at 
which the mean square of the density contrast at time t is unity, {S'^){Rni{t),t) = 1.) 
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If the universe followed the Einstein-de Sitter behaviour, the evolution would enter the 
saturation regime around 10-100 billion years. 

Once structure formation has started, part of the universe is in a constant state 
of collapse, and part is always becoming more empty. Based on the Buchert equations 
(fT0l) -( fT3i) and the analysis of the two-region toy model, we would expect the backreaction 
in the real universe to be strongest when the collapsing objects contract from having 
occupied the maximum fraction of volume. Neglecting evolution in number density, 
this naively occurs after the objects have reached their maximum size. The issue is not 
entirely clear, however, and it could be that the important factor is the slowing down 
of the growth of the collapsing regions as their size relative to the horizon becomes 
practically saturated. At any rate, regarding the coincidence problem, it is encouraging 
that these simple arguments lead to a time which is in rough agreement with the era 
when acceleration has been observed. 

There is another reason to think that acceleration will not occur with the formation 
of the first bound objects, but later in the history of the universe. If we set up the 
two-region toy model discussed in section 13.11 so that the initial volumes and absolute 
magnitudes of the density perturbations in the overdense and underdense region are 
equal (which is arguably closer to the real situation at early times), the expansion 
goes smoothly from the Einstein-de Sitter behaviour (a oc t^^^) to the empty universe 
behaviour (a oc t). There is no slowdown period and no acceleration, as the variance 
of the expansion rate is too small. The calculation that was done in the toy model is 
more representative of structure formation in the late-time universe, when voids already 
occupy a large fraction of the volume. 

So, it is plausible that acceleration only occurs once structure formation takes place 
in an environment dominated by voids. In view of this it is, again, encouraging that 
simulations and analytical work show that voids grow to fill practically the whole space, 
with the void distribution dominated by voids of a characteristic size, which grows as 
larger structures become non-linear and smaller voids merge [93, 119, 130, 131]. Not all 
analyses of simulations agree on the distribution of void sizes, evolution of the fraction 
of space occupied by voids or the fraction of space in voids today. In [94] the growth of 
the void volume fraction is clear, but the value today is only ^ 0.6, whereas the authors 
of [132] find that the fraction reaches ^ 1 at present day. Some of the simulations use 
the ACDM model, where the present day is singled out by construction, so care must be 
taken in applying the results to a dust-dominated universe. Nevertheless, given that the 
contribution of the cosmological constant in small at high redshifts, one would expect 
the feature that voids do not dominate early on to remain valid, while the behaviour at 
low redshifts should be checked more carefully. 

These issues should be looked at in detail both analytically and with simulations. 
However, while simulations are a useful complement to observations and analytical work 
in determining what the structures present in the universe are actually like, finding out 
about backreaction with simulations is not straightforward. Because simulations use 
Newtonian gravity and periodic boundary conditions, backreaction vanishes identically. 
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as discussed in section 12.21 Thus the effect of inhomogeneity and anisotropy on 
the mean expansion rate would not be visible in simulations even if the background 
evolution were dynamically adjusted to take into account the structures, instead of 
being predetermined. For work on the local effects of inhomogeneities and anisotropies 
on the expansion in a backreaction context, see [42]. 

The backreaction conjecture. Let us summarise the physical picture of backreaction- 
driven acceleration. As perturbations become non-linear, overdense regions collapse and 
form stable structures, and underdense regions form voids which become ever emptier. 
The geometry of the hypersurfaces of constant proper time is no longer everywhere 
perturbatively near homogeneous and isotropic, but a fraction of space is taken up by 
non-linear structures in various stages of expansion or collapse. A given proper time no 
longer corresponds to a single expansion rate, but to a distribution of expansion rates. 
Nevertheless, the description of the expansion in terms of an overall scale factor remains 
valid, because the individidual non-linear regions are small compared to the scales over 
which we are measuring observables. This distribution evolves until it reaches a self- 
similar regime where the regions grow at the same rate as the visual horizon. This 
happens when the universe is around 10-100 billion years old. 

As light rays coming to us pass through hypersurfaces of proper time, the fraction 
of space in each stage of expansion (or collapse) that they encounter changes as new 
non-linear regions emerge and old regions evolve. The conjecture is that the fraction 
of space in collapsing objects first grows, diminishing the average expansion rate. As 
these regions collapse and their contribution to the expansion rate is overcome by that 
of underdense regions, the average expansion accelerates, as demonstrated with the toy 
model in section 13.11 (Note that this mechanism relies on the structures we know to 
exist in the universe, as opposed to speculation of a globally inhomogeneous state [70].) 

This backreaction conjecture should be quantitatively checked by studying realistic 
models which can be compared with observations. However, let us briefly discuss what 
we can say regarding compatibility with cosmological observations on the basis of the 
Buchert equations and the qualitative description we have outlined. 

3.3. Comparison with observations 

Spatial curvature and the density parameters. Acceleration driven by backreaction 
necessarily involves spatial curvature. It is instructive to look at the behaviour of 
the universe in terms of the density parameters. Dividing ( |TOl) and ( ITTi) by 3H^, we 
have [41,55] 



where fi„ = 87rGjv(p)/(3i/2), = -{'■'^'^ R) /{QH^) and fig = -Q/{QH^) are 



the density parameters of matter, spatial curvature and the backreaction variable. 
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respectively. As seen from the definition of Q in (1141) . the backreaction density parameter 
is just minus the relative variance of the expansion rate, plus the contribution of shear: 

According to a variety of observations, today we have 0.15 < Qmo ^ 0.35 [26]. A 
rough estimate of present-day acceleration from SNIa data is —1.2 < go ^ —0.3 [133], 
though the lower limit could be extended to at least —1.6, and the upper limit to above 
zero, so go could even be positive today [3, 14]. Strictly speaking, results from even 
fairly model-independent analyses which do not assume the FRW equations but only 
the FRW metric [3,4] cannot be directly applied to backreaction, because of the different 
behaviour of the spatial curvature. In general, the results of SNIa data analysis depend 
strongly on the parametrisation used [12], and the data seems to contain little model- 
independent information about the expansion beyond the fact that it has accelerated 
in the recent past [3,4]. The quoted value of go should be understood as a rough 
estimate which assumes that the universe has not decelerated strongly since it started 
accelerating, sufficient for our purposes. 

Given the values of Qmo^Qo: the relations fl2Tl) . fl22l) imply 0.9 < Qrq < 1.5, 
—0.7 < Qqo < —0.2. It might seem that such a large negative spatial curvature would be 
in clear conflict with CMB observations [129]. However, as with SNIae, CMB analysis is 
very model- and prior- dependent. In particular, the analysis leading to tight bounds on 
spatial curvature assumes that the equation of state of the negative-pressure medium 
which is driving the acceleration does not vary rapidly in time. In contrast, at the 
level of the average equations, backreaction looks like a medium which evolves from the 
dust equation of state to having negative pressure. Such behaviour allows a significant 
contribution from negative spatial curvature [13], as the evolution of the equation of 
state can cancel the geometrical effect of spatial curvature on the angular diameter 
distance. 

As with the SNIa data, the non-FRW behaviour of the spatial curvature means 
that the CMB analysis in FRW models cannot be applied to backreaction as is. Spatial 
curvature would be expected to be negligible before backreaction becomes important 
and then evolve rapidly, rather than going smoothly like a~^, and it could even change 
sign. This makes analysis of the CMB anisotropics more involved than in the FRW case, 
as the basis functions for slices of constant curvature cannot be straightforwardly used. 
However, if backreaction becomes important only at redshifts of order unity and below, 
one might expect the spatially fiat FRW framework to be a good first approximation, 
as the spatial curvature becomes significant only on the last leg of the journey of the 
CMB from the last scattering surface to us. At any rate, the issue of light propagation 
in such an inhomogeneous and anisotropic spacetime should be studied in the context 
of a realistic quantitative model. 

Dynamical spatial curvature is a qualitative feature which is unique to backreaction. 
It might be useful in distinguishing it observationally from FRW models, which can share 
many of the other features of backreaction, including slowdown before acceleration, 
oscillations between acceleration and deceleration and a low Hubble parameter (see 
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e.g. [134]). 

Large-scale variance. It has been argued that backreaction cannot explain the observed 
acceleration since the geometry of the universe is so smooth, i.e. so near the FRW metric 
[66,71]. However, after non-linear structure formation has started, the geometries of all 
regions of the universe are not perturbatively near each other, as measured by invariant 
quantities such as the scalar curvature. For example, the difference in the expansion 
rate between expanding and collapsing regions is of order one and their behaviour is 
qualitatively different; as we have noted, the associated breakdown of perturbation 
theory is well known in the context of the spherical collapse model. A perhaps even 
simpler example is that the space inside stabilised structures such as galaxies does not 
expand (or collapse), so the relative difference between the local expansion rate and the 
mean expansion rate is exactly unity. This difference between the local static metric and 
the global expanding metric is also well known, see [135] for discussion and references, 
and has recently been studied in the context of spacetime-dependent couplings [136]. 

It is true that despite the presence of non-linear regions, the density field is smooth 
when averaged over large scales, as positive and negative density perturbations cancel 
due to conservation of mass. For this reason the energy density of dust necessarily 
goes like (p) oc a~^, as shown by ( |T2l) . and can never lead to acceleration. However, 
as we have discussed, it is not the variance of the density but the variance of the 
expansion rate which contributes to the acceleration. As estimated above, a variance 
of ((6'^) — (6Y)/(6)'^ ^ 0.2 — 0.7 (plus the contribution of shear) today could explain 
the observations. The important question is then what is the fraction of space occupied 
by regions with non-linear perturbations. If it is small, the variance can be negligible, 
whereas if it is of order one, the expansion will accelerate, unless the effect of the large 
variance is overcome by the effect of shear. 

Analysis of voids in the two-degree Field Galaxy Redshift Survey found that 40% of 
the survey volume is taken up by voids [92], which were required to be very underdense, 
5 < —0.9, and to have a minimum radius of 10 /i~^Mpc. The mean density contrast 
of voids was found to be 5 = —0.94 ± 0.02 and the mean radius was 14.89 ± 2.65 
/i~^Mpc. These numbers are in rough agreement with analysis of earlier surveys, some 
of which found larger mean sizes and a bigger fraction of the volume in voids [90,91]. 
The estimate of the volume fraction is conservative, both because of the limit on 5 
and radius, and because the void-finding algorithm looks for spherical voids (some of 
which are then combined). Simulations have indicated that voids can have complicated 
shapes, and that large voids in particular are typically non-spherical [93, 137] (though 
see also [119, 131]). The minimum size may also be important. Analysis of voids in the 
Millennium simulation (using the ACDM model) using the same void-finding algorithm, 
but with a minimum radius of 6 /i~^Mpc, found the mean radius to be 10.45 /i~^Mpc, 
while the mean density contrast was as low as before, 5 = —0.92 [138]. The fact that 
the mean size is not only smaller than in [91,92], but near the minimum radius of those 
studies suggests that the contribution of small voids may be important. The prevalence 
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of small voids is also suggested by some simulations [94], though not by others [93]. In 
any case, one should be careful when comparing observations of the real universe with 
simulations which use the ACDM model. It also bears emphasising that there is no 
generally agreed definition of a void [137], so it is not straightforward to compare void 
properties between different studies. 

We can make a rough order of magnitude estimate for the lower limit of the variance 
implied by the density distribution found in [92] by assuming that the volume taken up 
by voids is a single homogeneous and isotropic region and the rest of the universe is 
another homogeneous and isotropic region. This is like the toy model in section 13. 
but with region 1 being a general underdense FRW universe instead of being completely 
empty. Putting in the numbers 6i = —0.94, vi = 0.4 and demanding that the mean 
density equals that of a spatially flat background gives {0"^) / {0)"^ ~ 1.08, or Qqq ^ —0.08. 
Taking, more consistently, the background to be the sum of the two regions, we get 
{6'^)/{6Y 1.04, or VLqq ^ —0.04. (In the latter case, we have to take = 0.45, as 
the toy model does not permit a lower value for the given density distribution.) These 
are not far from the the value required by observations, VLqq < —0.2. These are likely 
to be severe underestimates of the variance, as the essential collapsing regions have 
been smoothed over. A realistic estimate should also account for the shear, which could 
counter the effect of the variance. 

We can also estimate the evolution of the expansion rate on the basis of observations 
and ACDM simulations of Lyman-a absorbers at redshifts between 2 and 4.5 [140]. 
It seems that the Lyman-a absorbers follow the overall Hubble flow, and allow a 
fairly direct observation of the expansion rate. The distribution of expansion rates 
analysed in [140] displays the qualitative features which we have identified as crucial for 
acceleration. Most of the volume is underdense and expanding faster than the mean, and 
there is a small fraction of space which is collapsing. The fraction of volume taken up by 
both the rapidly expanding and the collapsing regions grows with decreasing redshift, 
and space becomes increasingly dominated by voids which expand faster than average. 
According to the backreaction conjecture, the volume fraction occupied by collapsing 
regions should decrease at low redshifts so that the mean expansion accelerates. As 
the redshifts analysed in [140] do not go below 2, it is not possible to check the 
behaviour in the redshift range where the SNIa data indicate acceleration. However, we 
do see from the simulations that the variance grows with decreasing redshift, implying 
increasing backreaction. At 2; = 3.8 we have (6'^)/(6')^ 1.05, at z = 3.4 we have 
{9^)/ {9)^ ^ 1.07, and at ^ = 2.0 we have {9'^) / {9^ ^ 1.26. In other words, Qq evolves 
from —0.05 to —0.26, plus the contribution of shear (which it has not yet been possible 
to observationally determine). 

Having voids dominate a large fraction of space not only contributes to the variance 
of the expansion rate and thus to acceleration, it is also necessary for boosting Ht above 
the Einstein-de Sitter value of 2/3. Fitting the parameters of the ACDM model to 
the three-year WMAP data gives Hq = 73.4+3^ km/s/Mpc and to = 13.73^0:^7 billion 
years [129], resulting in (neglecting correlation in the errors) Hoto = 1.03l;o;og. This value 
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is very model-dependent, and we can get a more robust estimate using the Hubble Space 
Telescope measurement of the Hubble parameter, which was recently revised downwards 
to Ho = 62.3 ± 5.2 km/s/Mpc [28] from the old value Hq = 72 ± 8 km/s/Mpc [139]. 
Taking = 13 ± 1 billion years, the new value of Hq gives 0.70 < i^o^o ^ 0.97. 

The updated lower value for Hq is easier for backreaction to accommodate, since 
Ht ^ 1 requires that almost all of the volume of the universe is taken up by voids which 
are almost completely empty. The volume domination of voids demonstrates what 
the decreasing spatial curvature which accompanies backreaction-driven acceleration 
implied by the Buchert equations means in physical terms. Having H^tQ < 1 is a 
natural outcome of the backreaction framework, related to the domination of the space 
by voids, whereas in the ACDM model it is a coincidence. Given that backreaction 
cannot raise Ht above 1 (assuming that matter can be treated as dust and vorticity can 
be neglected), one can rule out backreaction as an explanation for the acceleration by 
showing, in a model-independent manner, that Ht > 1. (In fact, it is enough to show 
that the local expansion rate satisfies the inequality 6t/3 > 1 somewhere.) 

The above estimates show that the qualitative picture of backreaction outlined 
on the basis of the Buchert equations and the two-region toy model is not in obvious 
disgreement with observations and simulations. In fact, naive quantitative estimates of 
the variance of the expansion rate give the right order of magnitude required for the 
backreaction acceleration mechanism to work. A serious comparison with observations 
and simulations will require both a realistic quantitative model for the backreaction 
as well as more careful interpretation of the observations (since the FRW metric 
cannot be used) and simulations (since in backreaction in the Newtonian limit vanishes 
identically as an artifact of periodic boundary conditions, and simulations with a ACDM 
background may be misleading). However, the naive estimates above show that the 
backreaction explanation for the acceleration is plausible given what we know about 
inhomogeneity and anisotropy in the universe. 

Deviations from the average. In a homogeneous and isotropic spacetime, it is 
guaranteed that a comoving observer will measure the average values for the expansion 
rate, energy density and other observables, since there is no difference between the 
average and local values. One would expect this to also hold for spacetimes which 
are perturbatively near FRW everywhere. However, in a universe which contains large 
inhomogeneities and anisotropics, one can ask what is the relevance of the averages for 
an observer making measurements at one particular location. This is essentially the 
question of why the approximation of looking only at an overall scale factor is valid, 
discussed as assumption 1 in section 12.11 There are two answers. The first is that most 
measurements of cosmological quantities are made indirectly via observations of the 
CMB, LSS or SNIa, which are mostly sensitive to large-scale properties of the universe, 
rather than local ones. The second, more pragmatic, argument is the proven success of 
the scale factor ansatz in fitting observations, and the fact that the size of non-linear 
structures is small compared to the visual horizon guarantees the consistency of the 
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ansatz. 

However, as the size of structures relative to the visual horizon is not entirely 
negligible and the variance of the expansion rate is large, one could expect to see 
directional differences in observables which are sensitive to the large-scale expansion rate 
and density distribution, such as the low multipoles of the CMB and the optical depth. 
One can think of this as analysing the statistical scatter around the average behaviour 
given by the Buchert equations. It is tempting to speculate that backreaction could 
thus link the observed acceleration to the directional large-angle anomalies seen in the 
CMB [5-9]. In the first-year WMAP temperature map, the southern hemisphere has 
an optical depth of r = 0.24l;Q Qy, while on the northern hemisphere the optical depth 
is consistent with zero, with an upper limit of r < 0.08 (in the frame which maximises 
the asymmetry) [6]. A better studied feature is the presence of a preferred direction 
in the low multipoles, which is correlated with the dipole and the ecliptic plane. The 
former suggests an effect related to the structures which are responsible for our proper 
motion with respect to the CMB [112] or systematics regarding the calibration with 
the dipole [8], while the latter points towards a systematical effect associated with the 
motion of the WMAP sateUite. 

From cosmology one cannot get a correlation with the ecliptic, and probably 
not with the dipole either. The ecliptic correlation could be a coincidence, as the 
low multipoles are still anomalous even if the ecliptic correlation is neglected, and 
its significance has gone down in the three- year WMAP data [7, 9] . Even neglecting 
the ecliptic correlation, it is not clear how backreaction would produce such a distinct 
preferred direction. On the other hand, backreaction could tie the directional anomalies 
with the observed lack of power at large scales, unlike most proposals, which explain 
the directional anomalies by adding a new contribution which aggravates the amplitude 
problem. Since the behaviour of the perturbations is directly related to the spatial 
curvature, which is in turn tied to the acceleration, a connection with directional 
variation of the Integrated Sachs-Wolfe effect (and the inferred optical depth) seems 
plausible. It is interesting that the template of the anisotropic Bianchi Vllh model 
fits the anomalies quite well [141]. The fitting template corresponds to a homogeneous 
but anisotropic universe with shear, vorticity and negative spatial curvature. Though 
the details of the Bianchi Vllh universe are different from what is expected from 
backreaction, this does show that a globally inhomogeneous and/or anisotropic universe 
can explain the anomalies. 

Whether any of the CMB directional anomalies are related to backreaction should 
be studied in a quantitative model of light propagation in an inhomogeneous and 
anisotropic spacetime. One would also expect to obtain a prediction for the directional 
variation of SNIa luminosity distances, which should be different from that of a 
linearly perturbed FRW model [23, 24], and which might be testable in the next 
generation of observations. Such directional variation could be another clear way, in 
addition to dynamical spatial curvature, of distinguishing between acceleration driven 
by backreaction and a homogeneous and isotropic medium with negative pressure (or 
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an equivalent modification of gravity). 
4. Conclusion 

Acceleration and inhomogeneity/anisotropy. The observational evidence for the 
acceleration of the universe is usually interpreted in the framework of linearly 
perturbed Friedmann-Robertson- Walker (FRW) models, which describe a universe that 
is everywhere almost homogeneous and isotropic. In the context of such models, a 
medium with negative pressure or modified gravity is needed to explain the observations. 
This leads to the coincidence problem: why has the exotic matter or strange gravity 
become important only recently? The most significant qualitative change in the universe 
around the era where acceleration has been observed is the formation of non-linear 
structures, so it seems a natural possibility that the observed deviation from the 
general relativistic prediction of the homogeneous and isotropic cosmological models 
with normal matter could be related to the breakdown of homogeneity and isotropy. 

The issue of cosmological homogeneity and isotropy has been extensively discussed 
over the years by George Ellis and collaborators, notably in the context of the 
observational program of cosmology [84-87]. One of the issues they have highlighted 
is that averaging and applying the field equations do not commute: in other words, 
the average properties of an inhomogeneous and/or anisotropic spacetime do not satisfy 
the Einstein equation. The task of finding the model that best describes the average 
behaviour of the inhomogeneous universe has been termed the fitting problem. 

The relativistic equations which describe the behaviour of average quantities in 
an inhomogeneous and/or anisotropic, but irrotational, ideal fluid universe have been 
derived by Thomas Buchert [41,47]. The Buchert equations show that it is possible for 
inhomogeneities and/or anisotropics to lead to accelerating average expansion in a dust 
universe, even though the local acceleration decelerates everywhere. They also show 
that the fraction of space occupied by non-linear regions is the determining quantity, 
not the size of the individual regions. Even when the average properties of space can be 
described in terms of an overall scale factor, the evolution of the scale factor does not 
necessarily follow the FRW equations. 

The possibility that inhomogeneities and/or anisotropics could lead to acceleration 
was studied in the context of linear perturbation theory in [54], and it was suggested 
that acceleration could be due to perturbations which have entered the non-linear regime 
but haven't yet stabilised. The possibility of acceleration via backreaction has been 
numerically verified [76,79,80]. However, the physics of how structure formation leads 
to acceleration and the question of why acceleration begins much later than structure 
formation have been unclear. 

We have now clarified these issues, which turn out to be intimately associated with 
the process of gravitational collapse. With a simple toy model, we have explicitly shown 
how overdense regions can first slow down the expansion, which then accelerates as these 
regions shrink and their contribution to the expansion rate decreases rapidly as they 
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collapse. 

We have also noted that the matter-radiation equality scale imprinted on the 
dark matter power spectrum leads to a preferred time for structure formation that is 
near the observed acceleration era. The typical size of collapsing structures relative 
to the visual horizon grows rapidly at the start of structure formation, but then 
slows down, saturating around 10-100 billion years. A naive look at observations and 
simulations of structure in the universe shows that the degree of inhomogeneity required 
for backreaction to yield acceleration is plausible. 

The backreaction conjecture. The backreaction conjecture for the acceleration is simple. 
According to the Buchert equations, large variance of the expansion rate leads to 
acceleration. The physical interpretation is simply that the relative volume taken up 
by the regions of space which are expanding faster will come to dominate over the 
slower expanding regions, so the average expansion rate will rise. Collapsing regions, 
i.e. regions with a negative expansion rate, give a large contribution to the variance, 
since they contribute positively to the mean square but negatively to the square of the 
mean. Structure formation involves gravitational collapse, and the size of the collapsing 
regions is largest at late times when acceleration has been observed. 

Such an explanation keeps the phenomenological successes of the FRW scale factor 
in fitting the observations, while avoiding the failure of the FRW equations, which 
has required the introduction of a medium with negative pressure or modified gravity. 
This is in contrast to models which propose explaining the observations by the effect 
of inhomogeneities on the propagation of light without having accelerated expansion 
[20,21], where the success of the ansatz that one needs only to look at an overall scale 
factor is accidental. 

In the context of FRW models, there have been attempts to connect the late- 
time acceleration to inflation (via making the same scalar field responsible for both), 
the era of matter-radiation equality (via a tracker field which reacts to the change in 
the background equation of state) and dark matter (via unified dark matter and dark 
energy). Backreaction involves a subtle link to all these issues. Inflation determines the 
initial amplitude of the density perturbations, matter-radiation equality starts the clock 
for structure formation, and the nature of dark matter determines the processed form 
of the power spectrum and the time of formation of the first generation of structures. 

With many previously unclear conceptual and qualitative issues settled, the task is 
now to build a realistic model and make quantitative estimates that can be compared 
with observations. The relevant aspects of observations and simulations should also be 
understood better. On the basis of general considerations we can already state that we 
should have Ht < 1, and that there should be observable amounts of spatial curvature 
(assuming that vorticity is negligible and that matter can be treated as dust) [65]. There 
may be a slowdown period preceding the acceleration, and the expansion may oscillate 
between deceleration and acceleration, but these issues have to be worked out in the 
context of a detailed model. 
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Note that there are no new fundamental parameters to adjust, and any unknowns 
are due to existing uncertainties about the power spectrum, the modeUing of structure 
formation and so on: the backreaction conjecture is eminently falsifiable. Backreaction 
analysis simply entails doing the usually implicit averaging in cosmology in a way that 
is both mathematically consistent and takes into account the structures that are known 
to be present in the universe, as has been advocated over the years in the context of the 
program of observational cosmology and related work. Backreaction offers an elegant 
possible explanation for late-time acceleration. Whether or not this possibility turns out 
to be realised, the effect of structure formation on the expansion rate should be carefully 
evaluated to solve the fitting problem and complete the program of determining the right 
equations for describing the overall behaviour of the universe. 
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